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We investigate the indirect interaction between two isolated impurities in a Luttinger liquid de- 
scribed by a microscopic lattice model. To treat the electron-electron interaction U the functional 
renormalization group method is used. For comparison we also study the U — case. We find that 
for a wide range of impurity parameters the impurity interaction V12 as a function of their separation 
r oscillates with decaying amplitude between being attractive and repulsive. For half-filling of the 
band and in a crossover regime between weak and strong impurities the interaction becomes purely 
attractive. For U — and independent of the impurity strength the amplitude of the interaction 
energy falls off as 1/r. For U > the decay for small separations and weak to intermediate impu- 
rities is governed by a U dependent exponent larger than —1, which crosses over to —1 for large 
r. The crossover scale depends on the impurity strength and U. We present simple pictures which 
explain our results in the limits of weak and strong impurities. We finally also consider attractive 
interactions U < 0. 
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I. INTRODUCTION 

The physical properties of interacting one-dimensional 
(ID) electron systems differ substantially from the 
generic Fermi liquid behavior of higher-dimensional sys- 
tems, offering a whole variety of interesting new effects. 1 
Replacing the Fermi liquid concept of quasi-particles, the 
Luttinger liquid (LL) phenomenology has proved to cap- 
ture the low-energy physics of a large class of models. 2 
The low lying excitations of the system are no longer de- 
scribed by quasi-particles, but rather by collective den- 
sity excitations. Accordingly, even the presence of a sin- 
gle isolated impurity can have strong effects on the prop- 
erties of the system<£ii£i£i£ 

For the case of two isolated impurities placed in 
a LL the linear conductance has been investigated in 
detail^ i 10 ' 11 ! 12 ! 13 Other aspects have not yet been dis- 
cussed to this extent. Here we study the indirect in- 
teraction between two impurities mediated by the elec- 
trons. A particular promising candidate for measuring 
such forces are impurities in fermionic "atomic quantum 
wires" realized with ultracold gases»« For a continuum 
model of noninteracting electrons with delta impurities 
it was shown that the impurity interaction V12 as a func- 
tion of the impurity separation r oscillates between being 
attractive and repulsive. Its magnitude decays as 1/r. 
The period of the oscillation is it/Icf, with fcp being the 
Fermi momentum^ The effect of a repulsive interaction 
was then taken into account using a field theoretical ef- 
fective low-energy model and considering the two limits 
of weak and strong bare impuritiesi 15 i 16 The weak impu- 
rity case was analyzed using linear response theory while 
for strong impurities methods developed in the context of 
quantum Brownian motion were applied^ In addition, 
for strong repulsive electron-electron interactions bound- 
ary conformal field theory was used to study the impurity 
interaction^ 

Here we supplement the earlier studies. We consider 



a microscopic lattice model, focus on small to interme- 
diate electron-electron interaction and apply a method 
which is nonperturbative in the strength of the impuri- 
ties. This allows us to study the crossover from the weak 
to the strong impurity behavior. For the two limits sim- 
ple pictures of the observed physics are presented. Some 
emphasis is put on the question whether the impurity in- 
teraction as a function of r continues to oscillate also in 
the presence of electron-electron interaction! 15 i 16 ' 18 We 
briefly discuss the case of attractive interactions. 

The method of choice for our calculations is the re- 
cently developed functional renormalization group (fRG), 
which has turned out to be a powerful tool in the study 
of low-dimensional interacting electron systemsj 19 i 20 We 
use an approximate truncation scheme that by compari- 
son to exact and numerical results was shown to provide 
a good approximation for small to intermediate electron- 
electron interaction. It was used for interacting systems 
of up to 10 7 lattice sitesJ^i 

The paper is organized as follows. In Sec. |TT] we in- 
troduce our model and in Sec. IIIII give a brief outline of 
the fRG procedure used. In Sees. IIVI and [V] we present 
our results for the r dependence of the impurity interac- 
tion considering noninteracting and interacting electrons, 
respectively. We establish contact to the former calcula- 
tions. Finally, in Sec. I VII our findings are summarized. 



II. THE MODEL 

As our microscopic model we use the lattice model 
of spinless fermions with nearest-neighbor hopping and 
nearest-neighbor interaction on a large but finite number 
N of lattice sites. To suppress the effect of the bound- 
aries the interacting chain is coupled to noninteracting 
semi- infinite leads via adiabatic contacts. They are real- 
ized by varying the interaction smoothly across the two 
contacts from zero to the bulk value U We are mainly 



interested in the half filled band case. 

The kinetic part of the Hamiltonian reads 
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where c] and c 7 - denote the fermionic creation and an- 
nihilation operators on lattice site j. We set t = 1 and 
use the hopping as the unit of energy. In addition, the 
lattice spacing is set to one. The electrons are assumed 
to interact on the bonds between the sites [l,iV] via a 
spatially dependent nearest-neighbor interaction 
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defining the interacting wire of interest. The operators 
on the right-hand side of Eq. |(3J) are the local density 
operators rij = acj shifted by —1/2 to assure that the 
interacting part of the system is half filled. To turn on 
the interaction smoothly we use the function 

arctan[s(j - j„)] - arctan[s(l - j s )\ 
u o,j+i - u — 7—7-7N Tvi -..^-_r-^ 7TT W 
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with s = j and j s = 56 for the left part [1, ^] of the 
interacting wire and a similar function for the right part. 
In combination with the semi-infinite leads this provides 
us with a wide region in the center of the interacting wire, 
which has a constant interaction strength U . Hardly any 
effects from the contacts can be detectedj 13 ' 21 

In the region of constant U we place the two impurities 
whose interaction we want to study. We mainly consider 
site impurities, modeled by the Hamiltonian (V a > 0) 



(4) 



and sketched in Fig. [T] but also use hopping impurities 
described by (0 < t a < 1) 



H.c. 



(5) 



For both types of impurities the impurity separation r is 
defined as r — ji — ji ■ 

The homogeneous model H = ffkin + flint with con- 
stant interaction between all bonds (that is also in the 
semi-infinite leads) can be solved exactly via a Bcthc 
ansata^ and exhibits LL behavior for all fillings and 
all U except for half-filling with \U\ > 2. 23 In the half 
filled case, the relation between the LL parameter K, 
which later on will become important, and the interac- 
tion strength U can be given in a closed for m 22 ' 23 
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FIG. 1: (Color online) Schematic plot of the lattice model 
with site impurities; in the case of hopping impurities t / 1 
for two bonds 



III. FUNCTIONAL RENORMALIZATION 
GROUP 

A. General setup 

The general formulation of the fRG for the generating 
functional T of the one-particle irreducible n-particle ver- 
tices starts by introducing an infrared cutoff A in the free 
propagator leading to a A dependent T A . We here use a 
sharp energy cutoff— By differentiating with respect to 
A one can derive an exact infinite hierarchy of coupled 
differential flow equations for the vertex functions i 19 i 20 

The adaption of this general scheme to our present 
problem, namely to a spinless inhomogeneous LL at tem- 
perature T = 0, and the approximations involved in de- 
riving a closed set of equations for the self-energy are 
described in great detail in Ref. |2l|. The self-energy is 
approximated as frequency independent T, A (iuj) — > S A . 
Defining G A (iuj) = [G (icj)-E A ]- 1 , with the Green func- 
tion Go obtained from the kinetic part Eq. (TTJ) of the 
Hamiltonian, the flow equations for the spatial (Wannier 
basis) matrix elements read 
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The 2-particle vertex was parametrized by a static 
nearest-neighbor interaction of strength U A j +1 which im- 
plies that the self-energy is a tridiagonal matrix. The 
flow equation © of the 2-particle vertex is the especially 
simple form for half-filling, derivable from the general 
formula at arbitrary filling^ The fRG flow leads from 
A = 00 down to A = 0, where the original system is 
recovered. At the end of the flow ^fji° present the fre- 
quency independent approximation for the self-energy. 
This approximation scheme to the full hierarchy of flow 
equations was successfully used to study various aspects 
of inhomogeneous LLsj 12 i 13 ' 21 

For our lattice model the matrix elements of the inverse 
of the full Green function at scale A, [G A ]~ 1 (iA) for j, j' S 
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[1, AT] are given by 

A ]7j,(<A)=iW 

i (iA - #M) i&^hi + - 4; 
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(10) 

The second to last term effectively accounts for the semi- 
infinite leads as an additional contribution to the self- 
energy^ 

The initial values at cutoff Ao ^> 1 of the fRG flow 
equations and |8|) in the case of site impurities are 
given by2i 



= E VaS jJa , 
a=l 

whereas for hopping impurities one finds 
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B. O-particle vertex and grand canonical potential 

Earlier papers using the fRG for ID systems mainly 
focused on the calculation of the self-energy and observ- 
ables which can directly be computed from it. In the 
present paper we are interested in the grand canonical 
potential (GCP) from which the energy of the indirect 
impurity interaction can be obtained. The connection 
between the GCP and the O-particle vertex (OPV) is es- 
tablished by replacing the normalizing grand canonical 
partition function 3f in the thermodynamical average by 
the grand canonical partition function of the free sys- 
tem, that is the system without interaction and impu- 
rities. This keeps the vertex functions of order greater 
or equal one unaltered and merely changes the physical 
meaning of the OPV as can be seen by looking at its def- 
inition in terms of the generating functional T: 24 
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with the external source fields <j> and (j>. Thus, the OPV 
describes the difference between the GCP of the full sys- 
tem and the one of the free system. 

For our sharp energy cutoff the flow equation of the 
OPV reads^^o 



1 
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were the trace and the matrix structure refer to the quan- 
tum numbers and we explicitly introduced the conver- 
gence factor e luJV , which is usually suppressed. The limit 



r] \ has to be taken at the end of the calculations. 
The validity of Eq. (flu]) is not restricted to the approx- 
imation scheme used here, which leads to a frequency 
independent self-energy. 

Since in a numerical solution the flow necessarily starts 
from a large but finite cutoff Ao, we have to take into ac- 
count the flow from A = oo to Ao in the initial values. For 
the self-energy and the 2-particle vertex this procedure is 
described in Ref. HJ. In this large u limit on can replace 
E A (iw) by its frequency independent part £ A [even if a 
more sophisticated truncation than the one leading to 
Eqs. 0-© is used]. Approximating Go by its high fre- 
quency behavior Go-jj' (ioj) ~ 5j, j' / (iu>) and expanding 
the logarithm we obtain with f2 A= °° = 

ftA,*! r dA sin ^ A ) T rS] A . (17) 
With the large A behavior of E A (see Ref. UK) 



where Ii t j-,k,l is the bare antisymmetrized interaction, one 
obtains performing the integrals over A and A' in the 
limit r] \, 



(19) 



Here W denotes all single-particle terms of the Hamilto- 
nian not taken into account in the free propagator Go<^ 
For our problem this yields the initial values 
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for site impurities and 
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(21) 
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for hopping impurities. After this step r\ can be set to 
zero in Eq. (fTB|) . Using the matrix identity det e A = e TrA 
Eq. (fl6|) simplifies to 
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IV. NONINTERACTING CASE 

A. Analytical results 

In the noninteracting case, we compute the energy of 
the indirect interaction of two site impurities by evaluat- 
ing the partition function. For finite temperatures T this 
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gives 

Jp n{l-G ;ij(^„)[^+^] 

+ V X V 2 [G 2 . jd (ioj n ) - G 2 . j>j+r (iu> n )] } , 

(23) 

where the product involves all Matsubara frequencies u> n 
and r € N denotes the distance between the two impuri- 
ties. At half-filling (with chemical potential fi — 0) the 



Green function of the noninteracting and impurity free 
system is given by 

G^{z) = -^(-\ + \yf^ . (24) 

The interaction energy between the impurities Vi 2 (r) 
as function of their separation r can be obtained as the 
difference of the GCP for infinitely separated impurities 
and the GCP for finite This leads to 
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where we took the limit T — > and replaced the sum 
over Matsubara frequencies by an integral. 

In the limits of weak V\ , V 2 "C 1 and strong V\ , V 2 ^> 1 
impurities and for sufficiently large r (that is for r 3> 
l/kp\ throughout this work we will not be interested in 
the behavior of Vi 2 at very small separations r < 1//cf) 
Eq. ([25]) with the Green function Eq. p4|) can further be 
evaluated. For small Vi, V 2 <C 1 we obtain 



V 12 = (-1) 
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and for large Vi,V 2 ^ 1 
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(27) 



In both limits the impurity interaction oscillates around 
zero with a periodicity of two lattice sites, which corre- 
sponds to -k /kp with the Fermi momentum kp = 7r/2 at 
half filling, and decays as 1/r with increasing impurity 
separation. For strong impurities the oscillation around 
zero is asymmetric, while it is symmetric for weak im- 
purities. Furthermore, increasing the impurity strength 
from weak to strong for fixed r the impurity interaction 
changes from being attractive to being repulsive (even 
separations) and vice versa (odd separations). 



B. Numerical results 

For weak and strong impurities the r-dependence of 
Vis is shown in Figs. [3] and [3J respectively. To obtain 
these results we inserted the Green function Eq. (f2"4")l into 
Eq. (|25p and performed the integral numerically (not re- 
lying on any approximations). For simplicity we chose 
Vi = V 2 . To cover a wide range of separations the data 
are shown on a logarithmic scale. We always plot two 



subsequent separations. In the figures the triangles indi- 
cate a separation of an even number of lattice sites, while 
the circles symbolize odd separations. 

We can define an effective exponent of the decay of 
Vi 2 as the logarithmic derivative of the difference of the 
amplitudes with respect to the separation. This exponent 
is shown in the insets of Figs. [2] and [21 Consistent with 
the asymptotic analytical results Eqs. |26)) and (|2T[) the 
interaction decays with the inverse of the separation. A 
1/r decay is also found for impurities of intermediate 
strength as exemplified in the inset of Fig. 0] obtained for 
V\ = 1 and V 2 = 4 (for a discussion of the main part of 
Fig. [H see below). 



C. Interpretation 

For weak impurities the 1/r decay and the symmet- 
ric oscillations with period ir/kp can be understood from 
the spatial dependence of the Friedel oscillations of the 
electron density induced by a single impurity — For dis- 
tances sufficiently larger then n/kF the change of the 
density oscillates symmetrically around zero with period 
Tr/kp and in ID dies off as the inverse distance from the 
impurity. In lowest order perturbation theory in the im- 
purity strength this leads to a potential of similar shape 
at the position of the second impurity and thus to the 
observed behavior of Vis ■ This argument can also be ap- 
plied to hopping impurities and we again [as in Eq. (I26|) ] 
obtain a symmetrically oscillating (attractive for even r, 
repulsive for odd r) Vis which decays as 1/r. 

Also for strong impurities a simple picture of the r- 
dependence of Vi 2 can be given. To this end we consider 
a simplified model of a half filled tight-binding chain with 
open boundary conditions which follows in the limit of a 
single infinitely strong impurity (site or hopping). A sec- 
ond infinitely strong impurity is placed at distance r from 
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FIG. 2: (Color online) Oscillations of the impurity interac- 
tion as function of the separation for V\ = V% = 0.1 in the 
noninteracting case. Triangles: even separations; circles: odd 
separations. Inset: exponent of the decay. 
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FIG. 3: (Color online) The same as in Fig. [2] but for Vi 
V 2 = 30. 
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FIG. 4: (Color online) Disappearance of oscillations for Vi — 1 
and V2 = 4 in the case of noninteracting electrons. Triangles: 
even separations; circles: odd separations. Inset: exponent of 
the decay. 



the first one. The ground state energy of this setup can 
easily be computed analytically leading to Eq. (|27p if the 
infinitely strong impurities are modeled as site impuri- 
ties, while it gives 



Vv. 



-T7T- for r even 

12r 



or 
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(28) 



for hopping impurities. We thus find that for this type of 
impurities the interaction V12 is attractive (repulsive) for 
even (odd) separations in the limit of weak as well as of 
strong impurities. 27 Simple numerics shows that the sign 
also does not change for hopping impurities of interme- 



diate strength. 

This analysis reveals that the 1/r decay as well as 
the prefactors 7r/6 and —n/12 follow from an even sim- 
pler consideration. V12 can directly be extracted from 
the 1/N correction AE of the groundstate energy of a 
tight-binding chain with open boundary conditions. A 
simple analytic calculation shows that AE is given by 
n/(6N) for odd N and -n/(12N) for even N. For even 
N the chemical potential /1 = lies between the last oc- 
cupied and the first unoccupied state. In this case the 
results from (boundary) conformal field theory can be 
used. Within this approach one obtains for a model with 
conformal charge c = 1 



AE 



2AN 



(29) 



where v c denotes the charge velocity^ This result also 
holds for U 7^ 0. For the tight-binding model at U = 
one finds v c = 2 and Eq. (JUJ) leads to AE = -n/(12N) 
obtained above by direct calculation. To deduce the im- 
purity interaction from these considerations for hopping 
impurities the chain length N must be substituted by r. 
For site impurities a length N corresponds to an impu- 
rity separation r = N + 1. This explains the interchange 
of "even" and "odd" in Eqs. {27]) and l[28"]). The reason 
for the 1 jr decay of V\i for strong impurities is thus the 
appearance of 1/N corrections to the groundstate energy 
of a finite system with open boundary conditions. 

Our results for V12 in the lattice model are equivalent 
to the ones found in the noninteracting continuum model 
with 5 impurities investigated in Refs.[l5| andllrjl. There 
a different interpretation is given in the limit of strong 
impurities. 



D. Other fillings and fine tuned parameters 

We here assumed the system to be half filled. Relaxing 
this does not affect the general behavior of the impurity 
interaction, namely the oscillation of the interaction en- 
ergy decaying with the inverse separation. Merely the 
periodicity is altered. The behavior in the two limits 
of weak and strong impurities can again be understood 
within the two simple pictures presented above. 

We next describe a behavior which we exclusively ob- 
served for site impurities and at half-filling. In a narrow 
crossover regime between the limits of strong and weak 
impurities the interaction becomes entirely attractive, 
but still oscillates around a decaying average value with 
amplitudes which scale as 1/r (not shown here). A pe- 
culiar behavior is found under the constraint V2 = 4/Vi. 
For such fine tuned impurity parameters the oscillations 
completely disappear. The latter is depicted in the main 
part of Fig. |U Note that the amplitude still decays as 
1/r as shown in the inset. 
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V. INTERACTING CASE 

To compute the interaction energy V12 between the 
impurities for U ^ 0, we numerically solve the flow equa- 
tions ([I])-© and ([22)) with the appropriate initial val- 
ues. To isolate the energy of the impurity interaction, we 
evaluate the GCP for different configurations, namely for 
V\ = V2 — yielding the extensive part of the electronic 
interaction, V\ = and V2 ^ at site of interest and vice 
versa, finally the GCP of the full system with V\ ^ and 
V2 7^ 0. The interaction energy V12 is then given by 



V12 = (fi — Ovi=v2=o) — (£V 1= o#v 2 
— (£V 2 =o#Vi — £Vi=v 2 =o) • 



n 



Vi=V2: 



=0) 



(30) 



For the case of hopping impurities the same equation 
holds, if one replaces V a by t a and keeps in mind that 
the impurity free case corresponds to t± = t% = 1. 

For U ^ we examine systems of 2.5 x 10 5 sites, placing 
the impurities in the interval [22500, 200000] well apart 
from the sections of the lattice in which the interaction 
is turned on and off. This enables us to calculate the 
interaction energy for impurity separations of up to 10 5 
lattice sites. We restrict ourselves to half-filling and first 
consider bulk interactions U <G (0,2], which correspond 
to LL parameters K e [0.5, 1) [see Eq. (JBJ)]. 



A. Repulsive electron-electron interaction 
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FIG. 5: (Color online) Oscillations of the interaction energy 
as function of the separation for Vi = Vi = 0.1 and inter- 
action U — 0.5. Triangles: even separations; circles: odd 
separations. Inset: exponent of the decay. 



0.01 



-0.01 



h -0.7 



5 



-1.3 



10 
* » • 



separation [sites] 

10 2 10 3 10 4 



10 



10 



10- 



10 10 
separation [sites] 



10" 



to 



FIG. 6: (Color online) The same as in Fig. but for Vi = 
V 2 = 10. 



We focus on equal V\ and V2 but again verified that de- 
viating from this restriction does not qualitatively change 
our main conclusions. Further down we briefly comment 
on the behavior for fine tuned impurity parameters as 
discussed at the end of the last section for U = 0. Then 
it will also become important to consider V\ 7^ V2. 

Figures [5] and [6] show the r dependence of the interac- 
tion energy for weak and strong impurities. We find that 
the oscillation of the interaction with a periodicity of two 
lattice sites (corresponding to half-filling) is robust even 
in the presence of the electron-electron interactionj 15 i 16 i 18 
As for U = the impurity interaction for weak site im- 
purities is repulsive for odd r while it is attractive for 
even r. For strong impurities the opposite holds. Apart 
from the oscillation the interaction energy decreases as a 
function of r and we extract the "effective" r-dependent 
exponent as in the last subsection. It is depicted in the 
insets of Figs. [5] and [6j For strong impurities a clear 1/r 
decay is found. For weak to intermediate impurities as 
in Fig. [5] the exponent is different from —1. It starts well 
above —1 for small r and tends towards it with increasing 
separation. 

For weak impurities V\ and V2 and separations r suffi- 
ciently smaller then a crossover scale r c (see below), the 
indirect impurity interaction can again be understood us- 
ing linear response theory. In this parameter regime the 
Friedel oscillations in the electron density induced by a 
single impurity behave as (—l) r /r 2K ~ 1 i 21 ^ 29 This leads 



to a potential of the same spatial dependence at the po- 
sition of the second impurity^ This argument also holds 
for hopping impurities. We thus expect that the oscilla- 
tions of V12 are symmetric and that the effective exponent 
of the decay is 1 — 2A, which corresponds to —0.72 for 
U = 0.5 [see Eq. |(B}]. Our results are consistent with this 
argument (see Figs. [S] and [71 in the latter figure U = 1 
leading to 1 — 2 A' = —0.5). The small deviation of the 
exponent at small r can be explained by our truncation 
of the exact hierarchy of fRG flow equations. This leads 
to an approximate AfRQ which is slightly smaller than 
the exact 

The linear response analysis cannot be used at large 
Vi, Vi and in analogy to the U = considerations the 
appearance of the 1/r decay as well as the oscillatory be- 
havior can be explained in terms of the finite size scaling 
of the ground state energy of a system with open bound- 
ary conditions. In lowest order perturbation theory in 
U, AE is given by -tt(1 + U/tt)/{\2N) for even N. This 
is consistent with the (boundary) conformal field theory 
result Eq. ([29]) as v c (U) = 2(1 + U/ir) to lowest order 
in U . For odd N we did not succeed finding a closed 
form expression for the leading U, 1/N correction of the 
ground state energy. It is easy to show numerically that 
AE > for odd N. As in the noninteracting case this ex- 
plains the repulsive (attractive) impurity interactions for 
site impurities with even (odd) separation and attractive 
(repulsive) V12 for hopping impurities with odd (even) 
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r. We note in passing that our numerics (lowest order 
perturabtion theory as well as fRG) shows that for odd 
N, the U dependence of AE cannot solely be expressed 
in terms of v c (U), as it is the case for even N. 

For U > linear response theory in addition breaks 
down for r 3> r c even if weak bare impurities are con- 
sidered. This is related to the fact that in a LL a single 
impurity is a relevant perturbation in the renormaliza- 
tion group sense. Thus, for distances sufficiently larger 
then r c , that is in the low energy limit, even a single weak 
impurity acts as a strong perturbation (flow towards the 
"cut chain" fixed point).— Consistent with this pic- 
ture for U > and increasing r we observe a tendency 
towards the exponent —1 in the decay of V12 for weak 
to intermediate impurities (as in Figs. [5] and [7]). Al- 
though the renormalization group flow of the impurity 
determines the scaling exponents and thus leads to the 
asymptotic l/r decay of V12, it does not affect the sign 
of the impurity interaction. The latter is fixed by the bare 
impurity strength (see Figs. O and [S]) and was discussed 
above. This observation shows that even for r 3> r c a 
system with weak to intermediate bare impurities is not 
completely equivalent to a system with strong bare im- 
purities. 

In Fig. [7] we show the effective exponent of the decay 
for interaction strength U — 1 and various values of the 
impurity strength V\ = V% = V. The complete crossover 
from the weak to the strong impurity behavior of the 
decay exponent cannot be covered for a single V . This 
follows from the V and K (that is U) dependence of the 
single-impurity crossover scaled 



r c oc 



(31) 



For V <C 1, r c grows exceedingly, even beyond the large 
system sizes we can treat using the truncated fRG. The 
K dependence of r c can be read off from Fig. [5J where 
the effective exponent for V\ = V2 = 0.1 and various 
values of the interaction strength is shown. The inter- 
action ranges from U = 0.5 to U = 2, corresponding 
to LL parameters from K w 0.86 to K ss 0.5 [see Eq. 
©]. For small to intermediate U < 1 the effective expo- 
nent starts close to the linear response result 1 — 2K and 
tends towards the asymptotic value of —1 on a length 
scale decreasing with interaction strength. The asymp- 
totic exponent —1 is reached for separations beyond our 
system size. For larger interactions 1 < U < 2 already at 
the smallest separations considered here a clear deviation 
from the linear response exponent is observed, while the 
asymptotic exponent —1 is reached for r ~ 10 5 . 

Similar to the noninteracting case in a narrow regime 
of site impurity parameters V\ and V% the impurity inter- 
action changes from being attractive (repulsive) for even 
(odd) r to the opposite. In the crossover region V12 be- 
comes purely attractive and one can again fine tune V\ 
and V% such that V12 does not oscillate anymore. 
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FIG. 7: (Color online) Effective exponent of the amplitude of 
the indirect interaction for a collection of impurity strength 
Vi = V2 — V and interaction [7=1. The "noise" visible in 
the curve for V = 0.01 originates from the fact that V12 is 
very small for small V. 
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FIG. 8: (Color online) The same as in Fig. [7] but for fixed 
Vi = V 2 = 0.1 and different U. 
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FIG. 9: (Color online) Oscillations of the interaction energy as 
function of the separation for site impurities Vi = V% = 1 and 
attractive interaction U = —0.5. Triangles: even separations; 
circles: odd separations. Inset: exponent of the decay. The 
"noise" results from the small values of V12. 



B. Attractive electron-electron interaction 

Finally, we mention attractive interactions U < 
(that is K > 1) not considered so far. In this case 
a single impurity is irrelevant in the renormalization 
group sense^£&£ which suggests that regardless of the 
strength of the bare impurities for r»r c linear response 
theory can be used. We then anticipate that the im- 
purity interaction oscillates around zero and for r 3> r c 
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decays with the linear response exponent 1 — 2K. Our 
numerical results are consistent with this expectation as 
exemplified in Fig. [9] for site impurities with V\ = V2 = 1 
and U = -0.5, leading to K(U = -1) = 1.19 and 
1 — 2K = —1.38. The small deviation in the exponent 
again results from the approximate nature of our trun- 
cated fRG equations. A similar behavior in the decay 
exponent of the Friedel oscillations induced by a single 
impurity was discussed in Ref.[2l| (see Fig. 15). For U < 
accurate results for the exponent of the decay of V12 are 
difficult to obtain as the impurity interaction becomes 
exceedingly small. 

VI. SUMMARY 

To summarize, we examined the indirect interaction 
between two impurities placed in a LL. We started from 
a microscopic lattice model, using the fRG in the in- 
teracting case. We observed, that for weak as well as 
strong impurities the interaction oscillates between be- 
ing attractive and repulsive decaying algebraically with 
the separation of the impurities. In the case of non- 
interacting electrons the exponent of this decay is — 1 
independent of the impurity strength. For interacting 



electrons and weak impurity strength the powerlaw de- 
cay of V12 starts at a larger U dependent value 1 — IK 
and tends towards —1 for large separations on a scale 
r c depending on interaction and impurity strengths. For 
weak impurities this behavior can be understood from 
the oscillatory decay of Friedel oscillations of the electron 
density induced by an impurity using linear response the- 
ory. In the opposite limit the oscillatory 1/r decay can 
be traced back to the finite size scaling of the ground 
state energy of a chain with open boundary conditions. 
Our approach supplements analytic calculations obtained 
for large interactions^ and those performed in the lim- 
its of large and small impurities. 15 ! 16 It reveals the full 
crossover from weak to strong impurities. For attrac- 
tive interactions the impurity interaction oscillates and 
asymptotically (for large separations) decays with the 
linear response exponent 1 — 2K. 
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